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We develop an approximation scheme for the quantum mechanics of A" D0- 
branes at finite temperature in the 't Hooft large- A" limit. The entropy of the 
quantum mechanics calculated using this approximation agrees well with the 
Bekenstein-Hawking entropy of a ten-dimensional non-extremal black hole 
with 0-brane charge. This result is in accord with the duality conjectured 
by Itzhaki, Maldacena, Sonnenschein and Yankielowicz. Our approximation 
scheme provides a model for the density matrix which describes a black hole 
in the strongly-coupled quantum mechanics. 
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1 Introduction 



Properties of black holes in quantum theories of gravity have intrigued physi- 
cists for many years. In recent years much progress has been made in un- 
derstanding some of these properties using string theory. However because 
string theory is usually formulated perturbatively about a given spacetime 
background, it has been difficult to obtained a unified and complete descrip- 
tion of black hole physics. In the past few years nonperturbative formulations 
of string theory have been proposed in terms of large N gauge theories. These 
formulations give new hope for understanding some of these fundamental is- 
sues. 

In particular, Maldacena's conjecture gives a promising arena to answer 
some of these questions. Maldacena's conjecture |3J relates string theories 
in anti-de Sitter backgrounds to conformal field theories in a large- N limit. 
However there is no free lunch. In this framework, seemingly obvious ques- 
tions on the string theory side are hard to formulate on the gauge theory 
side. Moreover, in the range of parameters where semi-classical string theory 
may be used to construct black hole geometries, the dual gauge theory is 
strongly coupled. This makes understanding black hole states very difficult. 

One way of dealing with these difficulties is to perform Monte Carlo sim- 
ulations of the gauge theory at strong coupling [| . In the present work we 
pursue a different approach, in which we begin with an ansatz for the density 
matrix which describes the black hole in the continuum gauge theory. 

The version of Maldacena duality to be considered here relates black holes 
in ten dimensions with 0-brane charge to supersymmetric gauged SU(N) 
quantum mechanics with sixteen supercharges [ffl. The metric of the non- 
extremal black hole is 



ds = a' 



-h{u)dt 2 + h-\u)du 2 + cl/2{g f^ N)1/2 dn 2 s 



where c = 2 7 7r 9 / 2 r(7/2) and gyM is the Yang-Mills coupling constant. The 
horizon of the black hole is at U — Uq, which corresponds to a Hawking 
temperature 



5/2 

T = {9 2 y M N)- l/2 ( ) (2) 



U 



The dual quantum mechanics is to be taken at the same finite temperature. 
The black hole has a free energy, which arises from its Bekenstein-Hawking 
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entropy [|J. 

PF = -2.52 N 2 [ - 2 ~ 1/3 ) (3) 



rp \ I .Ml 



Duality predicts that the quantum mechanics should have the same free 
energy. The supergravity description is expected to be valid when the cur- 
vature and the dilaton are small near the black hole horizon. This regime 
corresponds to the 't Hooft large N limit of the quantum mechanics, when 
the dimensionless effective coupling gy M N/T 3 is large. 

In this report we describe a set of approximations that can be applied to 
the quantum mechanics in the regime of interest. Using these techniques we 
calculate the finite temperature partition function of the quantum mechanics. 
Over a certain range of temperature our results can be well fit by a power 
law, 

0F « const.- 2.0 ^(_L^y\ (4) 



This is in quite good agreement with the black hole prediction (|3p. We 
believe this is the first nontrivial direct test of a strong/weak coupling du- 
ality that does not rely on supersymmetric non-renormalization theorems 
or special properties of BPS states. Although in the present paper we are 
primarily interested in the thermodynamics of the quantum mechanics, our 
approximation scheme should also be useful for addressing questions about 
the spacetime structure of the black hole, perhaps along the lines of || 0. 

The basic idea is to treat the 0(N 2 ) degrees of freedom of the quantum 
mechanics as statistically independent, using a type of mean field approxi- 
mation. This assumption is motivated by the overall N 2 dependence of the 
free energy (||). The approximation involves constructing a trial action So 
from the full action S. All quantities can then be systematically computed 
as an expansion in powers of S — So- The parameters in the trial action 
are fixed by solving a truncated version of the Schwinger- Dyson equations 
of the quantum mechanics. This procedure can be viewed as re-summing an 
infinite number of Feynman diagrams. Since we are interested in large-A 
behavior, we will only resum planar diagrams. Thus, in our approximation, 
the overall N 2 factor in the free energy d) and the appearance of gy M only in 
the combination gy M N is guaranteed. The crucial test of the approximation 
is to obtain the correct dependence of the thermodynamics on the effective 
dimensionless coupling g 2 s = gy M N/T 3 . 

The sort of approximation that we are considering has several attractive 
features, which we regard as a priori reasons to believe that it captures some 
of the essential physics of the strongly-coupled quantum mechanics. 
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• As mentioned above, the approximation automatically respects 't Hooft 
large- N counting. 

• The approximation partially respects the symmetries of the problem. 
More precisely, it respects all symmetries which act linearly on the 
fields. Thus our trial action will have Af = 2 supersymmetry and 
5*0(2) x 50(7) rotational symmetry (out of the underlying M = 16 
supersymmetry and 50(9) rotational symmetry). 

• The approximation is non-perturbative in the Yang-Mills coupling con- 
stant, and self-consistently cures the infrared divergences present in 
conventional perturbation theory. 

Before continuing there is one other issue we would like to comment on. 
The partition function of the full quantum mechanics contains an infrared 
divergence from the regions in moduli space when the DO-branes are far 
apart. This leads to a divergent contribution to the entropy with an overall 
coefficient O(N). From the supergravity point of view, this corresponds to a 
thermal gas of gravitons. This divergence may be regulated by putting the 
system in a finite box. The black hole entropy which is 0(N 2 ) can easily be 
made to dominate over the O(N) contribution. Our mean field approxima- 
tion computes the 0(N 2 ) piece by design, and this infrared divergence does 
not make an appearance. 

2 Gaussian Approximation for O-brane Quan- 
tum Mechanics 

In this section we sketch the application of the Gaussian approximation H] 
to gauged SU(N) supersymmetric quantum mechanic with sixteen super- 
charges. Further details will appear in ||. 

One key requirement is that supersymmetry, although softly broken by 
the finite temperature of the black hole, should not be broken explicitly by the 
approximation. To avoid such explicit breaking we adopt an unconstrained 
superfield formulation, in which supersymmetry acts linearly on the fields. 
This ensures that we recover exact supersymmetry in the zero temperature 
limit. We will use an M = 2 superspace so that only an 50(2) x G 2 subgroup 
of the 50(9) R-symmetry is manifest. First let us recall the superspace and 
supermultiplets that we are going to use; for more details on notation see 0. 

With M = 2 supersymmetry we have an 50(2) R-symmetry, with spinor 
indices a, ft — 1,2 and vector indices i,j = 1,2. The 50(2)^ Dirac matrices 
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Y a/3 are real, symmetric, and traceless. Given two spinors ip a and Xa, besides 
the invariant ip a Xa, one can construct a second invariant which we denote 

Af = 2 superspace has coordinates (t, 6 a ) where 9 a is real. We denote 9 2 = 
6 a 6 a , and define the supercovariant derivative 

The simplest representation of supersymmetry is a real scalar superfield 

<5> = <P + t^J a + f6 2 (6) 

containing a scalar field 0, its superpartner ip a , and an auxiliary field /. The 
gauge superfield denoted by T a has an expansion in 'linear' components as 

r Q = Xa + A 6 a + X'^dp + de aP 6 p + 2e aP \p6 2 . (7) 

The fields X % are physical scalars, while \ a are their superpartners, d is an 
auxiliary boson, Xa are auxiliary fermions, and A is the 0+1 dimensional 
gauge field. 

To construct our action we introduce a collection of seven adjoint scalar 
multiplets $ a transforming in the 7 of Gi C 5*0(9). The M = 16 supersym- 
metric SU(N) gauged quantum mechanics action is 

S sym = 4~ I dtd 2 6Tr (~V a FiV a Fi - ^V a $ a V a $ a - l -f abc <b a [^ $ c ] j 

^YM J V 4 ^ / 

where jFj is the field strength constructed from the gauge multiplet, V a = 
D a + T ai and f a b c is a suitably normalized totally antisymmetric G^-invariant 
tensor. 

We impose the supersymmetric gauge condition D°T a = 0, which sets 
dA /dt = 0, d = and A = ^dx/dt. This is a convenient gauge fixing, as 
this gauge condition helps make the approximation compatible with Ward 
identities [[|. To the SYM action we must add the corresponding ghost action 
(but no gauge fixing term) 



ghost = 4~ I dtd 2 6 Tr (D a CV a C) 
9ym J 



9ym 

where the ghost superfield C is a complex scalar superfield with Grassmann 
statistics. 
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We are interested in the finite temperature behavior of the quantum me- 
chanics. As usual we compactify the Euclidean time coordinate on a circle 
of circumference /3, which is identified with the inverse temperature. The 
bosonic fields have integer mode expansions, while the fermions have half 
integer modes; for example we write 



Note that in Euclidean space the zero mode of the gauge field, which we 
denote A 00 , survives as a physical degree of freedom (fluctuations in A are 
eliminated by our gauge condition). 

Our strategy is to construct a trial action So, which we use as an approx- 
imation to the full action S = S sym + .Sghost- The free energy can then be 
calculated as an expansion in powers of S — S ; note that such an expansion is 
non-perturbative in the original Yang-Mills coupling g\ M . Throughout this 
paper we use the following expression for the free energy in this scheme. [] 



Here (3Fq is the free energy of the trial action, and (-)o denotes an expectation 
value computed using So- Also Sni refers to cubic terms in the original 
SYM plus ghost action, and the subscript C denotes a connected correlation 
function. It is straightforward, though tedious, to compute higher order 
terms in the expansion of j3F. In principle this could be used as a check on 
the validity of the approximation. 

We make the following ansatz for the trial action: 



Here all fields (except the gauge field) appear in Gaussian form. Indeed this 
is the most general Gaussian action which is quadratic in the fundamental 
fields. This means that the trial action can respect all symmetries which act 
linearly on the fundamental fields. 

The gauge field must be treated in a special way, owing to its periodicity 
properties. To do this we have introduced the timelike Wilson loop operator 
U, which can be expressed in terms of the gauge zero mode. 




PF « f3F + (S- S )o 




(8) 




jj — p e */c i A ° = e ^ A oo 



1 This quantity can be identified with the two-loop 2PI effective action of 
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This makes it manifest that at finite temperature Aq ~ Aq + 27[//3 is periodic. 
As a trial action for the gauge field we have adopted the unitary one plaquette 
model action. As A varies the trial action goes through a Gross- Witten phase 



transition at A = 2 1 10 



The key step in the approximation is to find a closed set of "gap" equa- 
tions for the dressed propagators of, Af, ... appearing in (|^). Again, the 
gauge field must be treated as a special case. All other propagators are 
obtained by demanding stationarity of the estimate ([8]) for the free energy. 
Up to contributions from the gauge field, it can be shown this procedure 
correctly re-sums all one-loop self-energy corrections to the propagators. 

The gap equation for A is obtained from the Schwinger-Dyson equation 
for (U) that arises from the change of variables U — > gU with g G SU(N). 
Demanding that this equation hold with respect to the one-plaquette measure 
yields 

This equation resums one-loop corrections to the Wilson loop, in the same 
sense that stationarizing d§|) resums one-loop corrections to the propagators. 
At large N the terms on the right hand side factorize into a gauge field 
correlator times matter field correlators; the terms involving the gauge fields 
may be computed using the results of [p!Cf| . 

As (§), fllOD are somewhat lengthy expressions we will not present them 
here. Rather we report on the solution to these equations in the next section. 



3 Numerical Results 



The gap equations can be solved numerically, using the methods discussed in 
Appendix B of ||. The basic strategy is to start at high temperature, where 
the gap equations can be solved semi-analytically, then use Newton-Raphson 
to solve the gap equations at a sequence of successively lower temperatures. 
As mentioned in the introduction, g\ M and N only appear in the combination 
gy M N/T 3 . Henceforth we choose units for T which effectively scales 9ym^ 
to 1. 

In principle the resulting Gaussian action contains a great deal of infor- 
mation about correlation functions in the quantum mechanics. But in this 
section we will just concentrate on the behavior of three basic quantities: the 
free energy, the Wilson loop, and the mean size of the state. 
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At high temperature, where the gauge theory is weakly coupled, we find 
that the free energy of the system is 

/3F = 61og/? + C>(l). (11) 

This result can be obtained analytically: the gap equations are dominated 
by the bosonic zero modes, and the free energy is dominated by f3F$. 

In general, for a weakly-coupled theory in + 1 dimensions, one would 
expect the free energy to behave like log/?. But note that, even though 
the gauge theory is weakly coupled at high temperature, the perturbation 
series is afflicted with IR divergences. Thus, to determine the coefficient 
of the logarithm (which depends on the value of the dynamically generated 
IR cutoff) one must re-sum part of the perturbation series. This is a well- 
known phenomenon in finite temperature field theory |12| . In any case, we 



expect a priori that the Gaussian approximation gives good results in the 
high temperature regime. 

As the temperature is lowered the behavior of the free energy changes: at 
(3 ~ 0.7 we find that it begins to roll over and fall off as a non-trivial power 
of the temperature. In the range 1 < (3 < 4 the numerical results for the free 
energy are well fit by 

(3F ~ —0.79 — 2.0 f3~ 1 ' 7 . (12) 

This fit to the numerical results is illustrated in Fig. 1. Note that super- 
symmetry is crucial in making such power-law behavior possible. Without 
supersymmetry the free energy would behave as (3F ~ (3Eq in the low tem- 
perature regime (3 > 1, where E is the ground state energy of the system. 
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We obtained ([12]) by performing a Levenberg-Marquardt nonlinear least- 
squares fit to 75 numerical calculations of the free energy, carried out in 
the temperature range 1 < /? < 4. To estimate the uncertainty in the best 
fit parameters we varied the window of (3 over which the fit was performed 
(fitting over the ranges 2 < (3 < 4 and 1 < (3 < 3), which leads to: —0.79 ± 
0.06, -2.0 ± 0.1 and -1.7 ± 0.2. 

It is quite remarkable that the power law (0) is in excellent agreement 
with the semiclassical black hole prediction 0J] 

f3F = -2.52 /T L8 ° . (13) 

The exponents differ by 6% while the coefficients of the power-law differ by 
26%. (An additive constant appears in the Gaussian approximation for the 
free energy. We will generally ignore this 'ground state degeneracy', since 
it seems to be an artifact of the Gaussian approximation when applied to 
systems with a continuous spectrum. Similar behavior was noted in 0.) 

As we go to still lower temperatures, we find that the energy 8{j3F)/df3 
calculated in the Gaussian approximation begins to drop below the energy 
of the black hole. In fact the Gaussian energy becomes negative around 
(3 = 5.8. Ultimately, as (3 —* oo, the Gaussian energy does asymptote to 
zero, as required by the M = 2 supersymmetry which is manifest in the 
approximation. But a negative energy clearly reflects some problem with the 
approximation. 

Fortunately, we can be rather precise about exactly where the approxima- 
tion is going wrong: the difficulty is with the Schwinger- Dyson gap equation 
we have been using to fix the value of the one-plaquette coupling A. Although 
we do not know how to write down a better gap equation for A, we can give 
a prescription for fixing A, that will allow us to obtain reasonable results at 
much lower values of the temperature. This may be regarded either as a 
check on our understanding of why the approximation is breaking down, or 
as a way of building a model for the black hole that can be used at lower 
temperatures. Our prescription for fixing A is simply that, when (3 > 2.5 
(the midpoint of our range 1 < (3 < 4), we choose A so that the free energy 
is given by (|12"D. 

The energy E = d(f3F)/d(3 calculated with this prescription is shown 
in Figs. 2 and 3. The behavior of A given by this prescription is shown in 
Fig. 4. Note that A increases monotonically with (3. A Gross- Witten phase 
transition takes place when A = 2; this value is reached at (3 = 7.8. Thus a 
phase transition takes place as the system moves into the supergravity regime 

i. 

By adopting the prescription of fitting j3F to a power law, we cannot say 
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Figure 2: Energy vs. (3. For f3 > 2.5 fixing A by fitting j3F to a power law 
leads to the solid middle line, while the Schwinger-Dyson gap equation for 
lambda leads to the lower dashed line. The upper dot-dashed line is the 
semiclassical energy of the black hole. 
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Figure 3: Energy vs. (3. Same as Fig. 2, but plotted on a log-log scale. 
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Figure 4: The one-plaquette coupling A vs. (3. The Gross- Witten transition 
occurs when A = 2. For (3 < 2.5 we use the Schwinger- Dyson gap equation 
to determine A. For (3 > 2.5 the Schwinger-Dyson gap equation gives the 
dashed line, while fitting (3F to a power law gives the solid line. 

anything about the order of the phase transition. If one takes the Schwinger- 
Dyson result for A seriously, then the Gross- Witten transition occurs at f3 — 
14.2, and is second order (the second derivative of the free energy drops by 
0.01 in crossing the transition). 

Our prescription for choosing A begins to break down around (3 = 14, 
as we find that A rapidly diverges as (3 approaches 14. By itself, this is 
not necessarily a problem: infinite A simply means that the Wilson loop is 
uniformly distributed over U(N). But unfortunately, we do not have a good 
prescription for continuing past this temperature. Evidently some of the 
other gap equations (not just the gap equation for A) start to break down 
at this point. Note that this breakdown does not occur until well into the 
strong coupling regime, as an inverse temperature f3 — 14 corresponds to an 
effective gauge coupling g^ s = (3 3 ks 3 x 10 3 . 

Finally, let us comment on the average 'size' of the state. In our approxi- 
mation the scalar fields X ! (r) and a (r) are Gaussian random matrices, and 
their eigenvalues obey a Wigner semi-circle distribution. We can define the 
size of the state in terms of the quantities 

Scalar = ^Tr < (X*(t)) 2 > , 

Cuge = ^Tr((r(r)) 2 ) . (14) 

The radius of the Wigner semi-circle, given by 2\/~R? ) is shown in Fig. 5. 
Note that the radius stays fairly constant in the region corresponding to the 
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Figure 5: Range of eigenvalues (radius of the Wigner semi-circle) vs. f3. The 
upper solid red curve is for the scalar fields in the scalar multiplets; the lower 
solid blue curve is for the scalar fields in the gauge multiplet. The dashed 
black curve is the Schwarzschild radius of the black hole. These results were 
calculated with (3F fit to a power law for (3 > 2.5. 



black hole. However, because the superfield formalism we are using does not 
respect the full £0(9) invariance, the radius measured in the scalar multiplet 
directions is not the same as the radius measured in the gauge multiplet 
directions. At (3 = 14 we find 

2-Rscalar = 1-81 2_R gauge = 0.80 . 

This shows that, as expected, the trial action does not respect the under- 
lying SO(9) invariance. Nonetheless, the trial action may provide a useful 
approximate description of the black hole density matrix in the supergravity 
regime. 

In Fig. 5 we have also plotted the Schwarzschild radius of the black 
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holef] Uo/2tt = 1.76 (3~ 2 ^. Note that, as the temperature decreases, the 
Schwarzschild radius becomes much smaller than the radius of the eigenvalue 
distributions. It seems appropriate to identify the radius of the eigenvalue dis- 
tributions with the size of the region U <C (gy M NY^ 3 in which 10-dimensional 
supergravity is valid {§]. The Gaussian approximation should provide a good 
laboratory for studying the way in which the horizon of the black hole can be 
detected in the dual quantum mechanics, perhaps along the lines suggested 



4 Discussion 



To summarize, we have presented an ansatz for a trial action which captures 
some of the behavior of 0-brane quantum mechanics at large N and strong 
coupling. The parameters appearing in the trial action are chosen according 
to a set of gap equations which resum an infinite set of planar diagrams. The 
approximation automatically respects 't Hooft large-A counting, and also 
partially respects the supersymmetries and R-symmetries of the quantum 
mechanics. 

Our main result is that we obtained a non-trivial power law for the free 
energy, which is in remarkably good agreement with the black hole prediction. 
The power law was obtained by fitting numerical results calculated in the 
range 1 < (3 < 4; we would like to emphasize that in this range of temperature 
the approximation involves no arbitrary or adjustable parameters. We also 
gave a prescription for fixing the expectation value of the Wilson loop, which 
allowed us to extend our results to (3 = 14. This is well into the strong 
coupling regime, where the supergravity description of the system should be 
valid. 

Although our prescription for fixing the Wilson loop cannot be regarded 
as entirely satisfactory, it does let us build a model for the black hole density 
matrix in the strongly-coupled quantum mechanics. This model could be a 
basis for further studies of black hole properties. It would be particularly 
interesting to address questions of spacetime locality and causality from the 
gauge theory point of view, perhaps along the lines suggested in || [7| . 
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